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ANALYSIS OF OSCILLATORY PRESSURE DATA 


INCLUDING DYNAMIC STALL EFFECTS 

Ry Franklin 0. Carta 
United Aircraft Research Laboratories 


SUMMARY 


The dynamic stall phenomenon was examined in detail by analyzing an 
existing set of unsteady pressure data obtained on an airfoil oscillating 
in pitch. Most of the data were for sinusoidal oscillations which penetrated 
the stall region in varying degrees, and here the effort was concentrated 
on the chordwise propagation of pressure waves associated with the dynamic 
stall. It was found that this phenomenon could be quantified in terms of 
a pressure wave velocity which is consistently much less than free-stream 
velocity, and which varies directly with frequency. It was also found that 
even when the stall region has been deeply penetrated and a substantial 
dynamic stall occurs during the downstroke, stall recovery near minimum 
incidence will occur, followed by a potential flow behavior up to stall 
inception. 

\ 

A small portion of the original experiment was performed with ramp 
motions in which linear variations of incidence angle with respect to time 
produced large regions of constant angular velocity. Here it was found 
that both stall inception and stall recovery were dependent on the 
instantaneous value of angular velocity for high incidence angles, but 
showed some dependency on past history for low incidence angles. 

A Fourier analysis was also performed on all data, and tabulations of 
pressure amplitude and phase angle harmonics are included. 



INTRODUCTION 


It has been known for many years that a torsionally flexible airfoil 
at high angle of attack is susceptible to a self-induced, single-degree-of- 
freedom instability commonly referred to as stall .flutter. As early as 193,8 , 
Bratt and Scruton (Ref. l) were conducting tests to measure moment , hysteresis 
in pitch and were using the concept of work per cycle around the moment loop 
to determine stability derivatives . Definitive and comprehensive studies 
were carried out later by Victory (Ref. 2) in Great Britain and by Halfman 
and his colleagues (Ref. 3) in this country. In these and in many other 
early experimental investigations .the measurements were necessarily confined, 
by the limitations of instrumentation then available, to gross measurements 
of unsteady lift, moment, and aerodynamic damping ratio for harmonically 
oscillating airfoils, and to gross measurements of response amplitude for. . 
freely fluttering airfoils. Hence, for many years the problem of stability 
at high incidence could only be treated phenomenologically, and the details 
of the unsteady flow over the airfoil could not be determined.. 

In 1957 Rainey published a report (Ref. 4) which described, among other 
things, the use of miniature pressure transducers distributed chordwise over... 
the airfoil to measure unsteady load distribution. Rainey's work is primarily 
concerned with the integrated unsteady forces and moments and his treatment, of 
the unsteady pressure measurements is limited to a brief discussion .in an . 
Appendix. Similar dynamic stall experiments were performed in the late . 

1950's (Ref. 5) in which NACA differential pressure transducers distributed , t 
along the chord were used to measure the unsteady pressures on a„.two-dimensional 
NACA 0012 airfoil oscillating in pitch. These data were integrated to yield 
unsteady lift and pitching moment which, in turn, were applied to various 
analyses of helicopter rotor blade dynamic stability as described in Refs. 6 
and 7* As in Ref. 4 the pressure data were primarily a means to obtain gross 
forces and moments, and only a cursory examination of unsteady chordwise 
loading was performed . A series of tests were then performed by Liiva et al 
(e.g.. Ref s . 8 and 9) in which unsteady pressures were again recorded and 
integrated to yield forces and moments . However, by this time the available 
recording and computing equipment had advanced to the point where it was 
feasible to also perform detailed analyses of the 0 basic pressure data and, 
for example, plots of pressure time histories are found in Ref. 8 and tabula- 
tions of harmonic pressure amplitudes and phase angles are found in Ref. 9* 

Within the past several years a serious effort has been made by many, 
investigators to explain the dynamic stall phenomenon in more detail, both 
analytically and experimentally. Of particular interest in these . studies has 
been the dynamic stall delay and the associated unsteady chordwise load, 
distribution on an airfoil moving through the steady-state stall regime with 
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positive angular velocity. In a recent theoretical study Ham (Ref. 10 ) was 
able to explain many of the features of stall delay and moment reversal on 
an oscillating airfoil operating above the stall angle by postulating a 
continuously shed vortex sheet from the airfoil leading edge. However, he 
was unable to derive a criterion for the delay between the time the steady- 
state stall angle is exceeded and the time the vortex sheet begins to be 
shed, nor was he able to predict reattachment of the separated flow (and 
hence closure of the moment loop) for decreasing angles of attack. In 
Ref. 11 Isogai indicated that the delay is associated with the formation and 
movement of an unsteady separation bubble. However, one remaining unknown 
in his investigation is the reason for the movement of the separation bubble 
against the adverse pressure gradient without bursting until well beyond the 
steady-state stall angle. These concepts are further discussed by McCroskey 
and Fisher in Ref. 12, who measured the chordwise motion of a shed vortex on 
a model rotor blade, and by Ham in Ref. 13, who speculates on bubble dynamics, 
bubble bursting, and dynamic stall as it relates to the adverse pressure 
gradient on a pitching airfoil. 

In a related analytical investigation (Ref. 14) it is shown that the 
unsteady pressure gradient over the forward portion of an oscillating airfoil 
is less unfavorable than the steady pressure gradient, which can contribute 
to the delay in stall for a dynamic process. McCroskey, in Ref. 15, points 
out that the theory of Ref. Ik agrees well with his experimental measurements 
for laminar separation, but he. also shows that the theory can only partially 
account for the dynamic stall delay. Thus, he concludes that although the 
inviscid theory of Ref. Ik accounts for the observed potential flow behavior 
of the oscillating airfoil, the viscous phenomenon of dynamic stall cannot 
be resolved by the same means . 

Crimi and Reeves (Ref. 16) have made the first important attempt to 
solve this problem of dynamic stall using viscous equations along the airfoil 
in which an adequate definition of the boundary layer is obtained. Included 
in their analysis are the so-called strong interactions between the external 
or free-stream flow field and the flow through the separated region. This 
analysis predicts a dynamic stall overshoot in both lift and moment, and the 
resulting plots of these response functions versus incidence angle display 
the characteristic loop closure observed in the many experiments described 
earlier. However, these predicted results do not agree well with experimental 
data, as indicated in Ref. 15, largely because this first effort is limited 
by the assumptions and restrictions of the formulation, such as the use of 
linear airfoil theory to predict the inviscid pressure distribution, an 
empirical transition model based on an experimental Reynolds number correla- 
tion, and use of an integral technique to model the separation region. An 
improved analysis is currently being performed under Contract NASI- 115 68 for 
the NASA Langley Research Center, which uses nonlinear unsteady aerodynamics. 
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a finite difference viscous flow analysis in the separation region, and 
an analytical turbulence - kinetic energy approach to predict transition. 

It is obvious from the foregoing paragraphs that, although a great deal 
of work has already been done, a significant amount still remains. Further- 
more, as more detailed analytical studies are attempted there is an increasing 
need for comparably more detailed experimental results , particularly in the 
dynamic behavior of the chordwlse load distribution. In a recent experimental 
program (Ref. 17) a substantial body of unsteady chordwise pressure data 
was recorded on FM magnetic tape. As in many other comparable experiments 
the emphasis in Ref. 17 was on the integration of data to yield unsteady lift 
and mcanent. In addition, the unsteady pressure data were carefully preserved 
and the detailed analysis of these data constitute the substance of this 
report. The original experiment was performed over a wide range of incidence 
angles, both above and below static stall, and for a number of oscillatory 
frequencies. Consequently, the character of the results obtained herein 
varies widely, from potential flow through a transition region up to a flow 
dominated by separation. Many of the results of this study can be explained 
and understood within the framework of the present state of the art; however, 
in view of the complexity of the dynamic stall phenomenon the interpretation 
of these results is necessarily incomplete, and it is hoped that additional 
studies may be performed in the future to answer these still outstanding 
questions . 
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P, P 
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T, T„ 


At l" At 2 
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t , t 
u D 
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V , V 

W W, 


a 


a 


M 


angular velocity parameter, ca/2V 
harmonic pressure amplitudes 
differential pressure coefficient 
chord, meters 

real and imaginary parts of Theodorsen circulation function 
frequency, Hz 

reduced frequency parameter, cu/2V 
average and initial wave number 
.harmonic index number 

pressure phase angle, positive when pressure leads motion, deg 

2 

upper and lower surface pressure, newtons/meter 
period, sec 
time, sec 

time intervals, sec 
dimensionless time, utj2v 

dimensionless upstroke and downstroke time for ramp cam motion 

free-stream velocity, meters/sec 

normalized average and initial wave velocity 

instantaneous angle of attack, deg 

mean incidence angle, deg 


vn 


angular amplitude, deg 


average and initial reduced wave velocity 

r 

density, kilogram/meter 

pressure phase angle, positive when pressure leads motion, 
radians 

dimensionless chord position 
dimensionless chordwise interval 
frequency, radians/sec 
derivative with respect to time, sec 



TEST PROGRAM 


The oscillating airfoil test program was performed for the U. S. Army- 
in 1971 - 1972 under Contract DAAJ02-71-C-0003, and the details of the model, 
test facility, data system, and test procedure are fully documented in Ref. 17* 
Hence-, ‘/the only aspects of the test program that will he discussed here are 
those that are germane to the analysis of the unsteady pressure data. In 
particular. Table II of Ref. 17 lists all the test conditions that were 
originally run, and the current study is ' confined to an -analysis of Items 
2, 7, and 8 from that table, which are reproduced here in somewhat greater 
detail in Table I. In the upper part of the table are the combinations of 
frequency and mean incidence angle at which test points were run for sinusoidal 
motion at an amplitude of a = 8 deg. The second-to-last column is the set of 
average frequencies generally used in this study, and the last column is the 
average reduced frequency, k = cu/2V, based on the average frequencies and on 
a blade chord of 0.127 m (5 in) and a free-stream velocity of II8.87 m/sec 
(390 ft/sec). The lower part of the table is for ramp motion, separated into 
forward and backward motions, denoted as such because the drive cam to produce 
the forward motion was reversed to produce the backward motion. This oscilla- 
tory schedule is completely documented in Ref . 17 and it is sufficient to note 
that over most of the motion the instantaneous angle of attack varied linearly 
with time (i.e., with constant angular velocity) over an incidence range of 
* 8 deg relative to the mean angle. For the forward ramp the upstroke velocity 
was twice that of the downstroke and for the backward ramp the upstroke velo- 
city was one-half that of the downstroke. Because this was a nonsinusoidal 
motion the frequencies listed here are the values of the fundamental (or 
one-per-rev) frequency of the drive cam. (Additional details of the ramp 
motions will be found in the last section of this report . ) 

Differential pressures were measured at ten chordwise locations, whose 
coordinates were chosen in accordance with a Gaussian integration procedure. 

In percent chord, these locations are 1.19; 6.15, 14.62, 25-83, 38*81, 52*38; 
65-36; 76-57; 85-04, and 90.00. In the original test program these data were 
recorded on FM magnetic tape and were then digitized and stored on digital 
magnetic tape by means of the data system described in Ref. 17- These 
digitized pressure time histories were the subject of the present analysis, 
as described below. 
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DATA REDUCTION 


Preliminary Data Preparation 


* j" 

‘X. r.J - 

- .i 
/ • •: 


As stated above , the unsteady pressure data for .each chordwise lQ ? cation 0 
were digitized and stored on digital magnetic tape. The digital sampling' 
rate was selected to yield a minimum of four • samples/cycle in the tenth 
harmonic of the fundamental frequency and it can be seen in Table III of 
Ref . IT that this minimum requirement was exceeded for most cases . This 
requirement of four or more samples/cycle in the tenth harmonic was specifi- 
cally chosen to permit a valid ten-term Fourier analysis to be made without ^ 
risk of higher harmonic distortion or encroachment into the lower harmonic . ... 

f *• 4" K i.' 'm- \ ,.'J ' * .’.i. -0 

data (e .g . , aliasing in which the higher harmonics are folded into the low 
frequency domain; see Ref. 18). . _ .. " . “ . 

Various aspects of the analysis were carried out on either the original 
time history or on an 8-cycle signal-averaged (or smoothed) time history 
(see Ref. IT and pp. 116-118 of Ref. 19). In each case the procedure used 
will be identified. At small values of mean incidence angle both methods 
gave virtually identical results. At large values of mean incidence angle, 
when stalling was present, the signal-averaging procedure had the effect of^" 
reducing the influence of random noise and enhancing the coherent or repeti- 
tive parts of the signal. 'A comparison of the raw data with signal-averaged 
data is shown in Fig. 1 for the 1.19 percent chord location at a.nqminal 
frequency of 98.5 Hz f or = 12 deg and 1 6 deg. It is seen that . although " 
there are minor differences between these plots, the essential characteristics 

1 r • ' • ■ T. ■?.<: ? ’ v- -. 

are retained in the smoothed plots. Further examples of the ability of the 
signal-averaged data to represent the raw data with satisfactory fidelity 
will be pointed out belpw. ......... 


Oil! 


nl::1 


, ; • Time History Displays . , ; . 

' ' . ' .* ' r * • **{*£? 

Time history displays were obtained by processing the raw time history . 
data described above on a PDP-6 computer and displaying the results on a 
Tektronics Model 6ll Storage Tube Oscilloscope (memo scope).. These displays w 
were then photographed for later analysis. " , 

The data were arranged within the computer such that at each instant of 
time a continuous function was generated to provide a polyn omi al data fit 
through the ten chordwise locations . Interpolation of the chordwise distri- 
bution was then possible at equally spaced chordwise intervals between the 
dimensionless chord positions X= 0.0119 and X= 0. 90 which were selected at " 

f \ • . v- Vi.'o : -.v: . f . > 

the convenience of the investigator (usually 33 intervals were chosen but any 
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integral number was possible). At each selected chordwise point a continuous 
time history was constructed with a sequence of straight lines between 
timewise points. Sufficient data were available to permit this type 
of fit and still yield a continuous behavior (from kO to 100 points per cycle 
as can be determined from Table II of Ref. 1?). A sample of this type of time 
history display is shown in the upper portion of Fig. 2 for = 3 deg and 
f = 98.5 Hz . Here the pressure difference coefficient, 

p (x,t) - P, (x,t) 

ACp(x,t) (1) 

Jpv 2 

has been plotted versus time for a number of interpolated chordwise stations, 
as discussed above. The upper curve, labled a( t), is the time history of the 
angle of attack, and the horizontal line, labled datum, is the zero datum for 
the leading edge pressure wave. The upper curve (below the a trace) is the 
pressure at the 1.19 percent chord location (denoted L.E.), the bottom curve 
is the pressure at the 90 percent chord location (denoted T.E.) and all of 
the curves between these two are the intervening plots at 32 (in this case) 
equally spaced interpolated chordwise locations . To avoid confusion, datum 
lines for chordwise pressures other than the leading edge pressure are not 
included . 

Because the vertical axis represents both pressure level and chordwise 
location the result is a pseudo-three-dimensional plot of the pressure time 
history surface from airfoil leading edge at the top to trailing edge at the 
bottom, with time increasing toward the right. The observer viewing this 
surface is positioned downstream of the trailing edge and above the zero 
pressure plane. 

The lower portion of Fig. 2 is a contour plot in which constant pressure 
levels have been outlined. Again the chordwise extent of the display is from 
the 1.19 percent chord to the 90 percent chord. This contour plot (and others 
to be shown subsequently) is best interpreted by comparing it with the time 
history directly above.it. (Note that in the contour plot the leading edge 
is at the bottom and the trailing edge is at the top.) The contour lines 
are crowded together in regions of large pressure gradient and are far apart 
in regions of small pressure gradient. Pressure levels are not explicitly 
included in the figures although the small superimposed numbers in the contour 
plots represent coded pressure levels, primarily intended for computer use. 

Further clarification of this combination of time history pseudo-surface 
and contour plot is shown in Fig. 3 which is a schematic version of the plots 
in Fig. 2. For simplicity, only a few of the contour lines have been retained 
in the bottom half of the figure and they have been sketched into the time 
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history surface in the top half of the figure. The appended numbers from 1 
to 5 relate the contour lines in the upper and lower portions of the figure 
to one another. Also, the pressure time history surface at the top of the 
figure for the first cycle has been truncated along contour line number 1 
to further relate the constant pressure surface in the upper portion of the 
figure with the constant level contour line in the lower portion. 

In the course of this study it was found that the pressure time histories 
were useful in visualizing the phenomenon of dynamic stall and the interplay 
between potential flow and stalled flow. The contour plots were useful in 
quantifying certain aspects of the chordwise propagation of pressure waves . 
This will be described at length below. All the pressure time histories and 
timewise contour plots described herein are for raw data only. 


Fourier Analysis 

The pressure difference time histories were also analyzed harmonically 
to yield the average value and the first ten Fourier harmonics (amplitude 
and phase angle) at each chordwise transducer location. Specifically, it 
was assumed that for sinusoidal motions the unsteady angular displacement 
and pressure difference coefficient time histories could be represented by 


«(t) = « M + SeltJt ( 2 ) 

ACp(x,t) = A q (x) + Y A n (x)e 1 ( ntJfc + (3) 

where X is a dimensionless chordwise coordinate normalized with respect to 
the airfoil chord, and where <£ n (x) is the phase angle in radians by which 
the nth harmonic leads the motion. For convenience a dimensionless time 
t* = ut/2Tr can be defined, and if P n (x) = 360 (j> n (X)/2iT is the phase angle 
in degrees, the real parts of Eqs . (2) and (3) become 


a(t*) = a + a cos 2wt* 
M 


(h) 


ACp (x,t*) = A q (x) + ^ A JX) cos 2 it 



( 5 ) • 
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Initially this ten-term harmonic analysis was performed on "both raw data 
(one cycle) and on signal- averaged data for all frequencies at a mean angle 
of attack of = 14 deg. It was found that for the first four harmonics the 
differences between the raw and signal-averaged results were small, and either 
of the two could he used to represent the original function. (For the higher 
harmonics the amplitudes are so small as to be within the error band of the 
data and the phase angles tend to become meaningless . ) A typical comparison 
between the two for = 14 deg and f = 50 Hz is shown in Figs . 4 and 5 for 
the amplitudes (note scale changes) and phase angles, respectively. At higher 
frequencies the two sets of data were virtually indistinguishable and at lower 
frequencies the agreement deteriorated only slightly. Therefore, because it 
was convenient to compute the harmonic coefficients for the time-averaged data, 
these have been used throughout the remainder of this study. Tabulations of 
all pressure harmonics will be found in Table II for the sinusoidal data and 
in Tables III and IV for the forward and backward ramp data, respectively. 
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ANALYSIS OF RESULTS 


Potential Flow Behavior 


The pressure time histories and timewise contour plots in Fig. 2 are 
typical of the potential flow behavior for low mean incidence angles and 
high frequencies . Here all the pressure waves are nearly sinusoidal, and it 
is seen that the trailing edge pressures lead the leading edge pressures. 
This is in agreement with potential flow theory. Specifically, the pressure 
phase angle, relative to the motion, can he derived from the formulas of 
Ref. 14 (with a suitable chordwise coordinate transformation and imposition 
of a quarter-chord pivot ) as 


<f>{x) = tan 


-1 


k 

kF + G - — ■ + 4kx' 

2 

2 2 

F - kG - 2k X 


(6) 


A comparison between the theory as predicted by Eq. (6) and the measured first 
harmonic phase angle distribution for a M = 3 deg and f = 31, 75 > and 98.5 Hz 
is shown in Fig. 6. The agreement is seen to be very good. 


For future reference it is important to note the behavior of the timewise 
contour plot in the lower portion of Fig. 2. It was stated earlier that the 
contours are crowded together in regions of high pressure gradient . It can 
also be seen that zero pressure gradient will cause any given contour line to 
attain its chordwise extremum. From another point of view, if this were a 
topographical map of a ridge line (in this case a pressure ridge) with a 
monotonically decreasing altitude, the constant altitude contour lines would 
form a nested family of curves having its axis aligned with the ridge . This 
is indeed the case in Fig. 2 in which the family of contour lines is seen to 
be inclined from upper left to lower right, indicating a trailing edge lead 
and a leading edge lag. 


Mixed Flow Behavior 

As shown in Fig. 18 of Ref. 17, static stall for this test program began 
at approximately 12 deg incidence (as measured from the departure from lin- 
earity) and was complete at approximately 16 deg. In the dynamic portion of 
the test program the oscillatory amplitude was a = 8 deg. Hence the maximum 
incidence angle during each cycle of motion for = 6 deg was 14 deg, for 9 
deg was 17 deg, etc. The unsteady results for = 6 deg were found to be 
nearly identical to those for o;^ = 3 deg, and thus were in good agreement with 
linear theory. This is not surprising because even at maximum incidence the 
static stall region has not been fully penetrated . 
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At = 9 deg the beginnings of dynamic stall are evident near the end 
of the upstroke of the motion. This is seen in Fig. 7 in which memo scope 
photographs of pressure time histories and pressure contours for = 9 deg 
and for f = 12.5, 31, and 98-5 Hz are reproduced. Before a reaches its 
maximum value the pressure time history attains a sharply peaked maximum 
value followed by a sudden drop in pressure. This sharp peak and sudden drop 
occurs first at the leading edge and propagates rearward along the chord, 
occurring at the trailing edge at some time later, and with less suddenness 
than at its inception at the leading edge . A deeper penetration into stall 
occurs at = 11 deg, as seen in Fig. 8 for f = 12.5, 50, and 98-5 Hz . 

Here the sudden drop in pressure is more precipitous and a depressed pressure 
level extends over a greater portion of the cycle, compared with the results 
in Fig . 7 for a M = 9 deg . 

In both of these figures there is evidence that a potential flow exists 
over the portion of the cycle from minimum a to a point just before the 
pressure peak occurs . This is best seen in the two right hand panels of 
Figs. 7 and 8 for f = 98-5 Hz . When these figures are compared with Fig. 2 
for = 3 deg it is seen that all three manifest the same behavior during 
the upstroke; viz., a nearly sinusoidal wave first appears at the trailing 
edge which appears to propagate toward the leading edge . An examination of 
the contour plots in all three figures (for f = 98-5 Hz) also shows a 
characteristic nesting of the family of curves with an axis that is inclined 
from upper left to lower right. Note, however, that the contour plots for 
ckm = 11 deg, f = 98.5 Hz in Fig. 8 exhibit a second set of extrema which 
are inclined from upper right to lower left, indicating a rearward propagation 
of pressure waves, in this case associated with dynamic stall. 

This mixed flow condition, consisting of potential flow during the 
upstroke and stall flow during the downs troke, is to be expected when the 
mean angle of attack is at or near the stall angle (such that excursions 
take place into both potential flow and stall flow). However, it is surprising 
to note the persistence of this mixed flow at mean angles of attack up to 
a M = 16 deg, as shown in Fig. 9 for f = 98*5 Hz. This condition appears to 
be frequency dependent and becomes much less pronounced at lower frequencies, 
but at high frequency the potential flow region can be observed in Fig. 9 as 
a sinusoidal wave in the upstroke portion of the time history, progressing 
from trailing edge to leading edge . It also appears in the portion of the 
nested family in the contour plots whose axis is inclined from upper left 
to lower right. Note that the contour plots in Fig. 9 show a diminishing 
potential behavior relative to the stall wave propagation as the mean angle 
of attack increases . Specifically, a prominent lobe to the left of the peak 
incidence angle which is associated with potential flow at = 12 deg is 
diminished at 14 deg and virtually disappears at 16 deg. 
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Stall Flow Behavior 


A series of memo scope photographs are presented in Figs . 10, 11, and 12 
to illustrate the effects of both frequency (over the entire range from 12.5 
to 98.5 Hz) and mean angle of attack = 12, 14, and 18 deg, respectively) 
on the pressure response during dynamic stall* In Fig. 10 for = 12 deg 
.the effect of increasing frequency appears to lessen the severity of the stall. 
At f = 12-5 H? a sudden stall occurs during the upstroke followed by recovery 
during the downstroke . In all succeeding parts of Fig. 10 recovery seems to 
be initiated at roughly the same part of the downstroke but the depth of the 
stall drop-off is diminished markedly as frequency increases . Similar behavior 
is observed in Fig. 11 for = 14 deg, but a higher frequency is needed before 
dynamic stall drop-off is reduced appreciably. This is seen by comparing the 
panels for f = 50 Hz in both figures . Although there is some effect of 
frequency at = 18 deg (Fig. 12) its influence is considerably less than at 
lower mean incidence angles, and the deep penetration into stall causes the 
same type of dynamic stall drop-off over the entire frequency range. A' 
comparison of the right hand panel of Fig. 9 for = 16 deg at f = 98*5 Hz 
shows a behavior more like = 14 deg than like = 18 deg. 

A recent paper by Martin and his co-workers (Ref. 20 ) is concerned with 
the dynamic stall phenomenon, and discusses, among other things, the rearward 
propagation of the dynamic stall cell along the chord of the oscillating 
airfoil. The description of Fig. 2 in Ref. 20 states that "the peaking of the ' 
leading edge pressure indicates that stall has begun". It is further indicated 
that the disturbance appears to progress toward the trailing edge, indicative 
of the passage of a vortex over the airfoil, and the appearance of another 
pressure disturbance is indicative of the formation of a secondary vortex. 

Multiple pressure waves are also observed in the present study, as seen 
in selected panels of Figs. 7 through 12. Here both secondary and tertiary 
pressure waves propagate rearward in the vicinity of the airfoil leading edge, 
although they appear either to merge with the primary wave or to attenuate to 
an undetectible size. It appears that in some instances (e.g., f = 50, 75 > 

98.5 Hz at Oty = 14 deg in Fig. ll), when a strong secondary wave is generated, 
the primary wave is retarded sufficiently (perhaps by some mutual interaction) 
to cause the two waves to coalesce, after which the rearward propagation of a 
single strong wave is observed. This obs-ervat ion will be of some value below 
in explaining the possible reasons for the appearance of a "kink" in the 
pressure ridge locus . 
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Wave Velocity Calculations from Chordwise Time Belay 

Measurements were made of the propagation rate of the pressure wave 
along the -chord lasing the contour plots prepared from the memo scope photo- 
graphs. Figure 13 depicts this schematically. The upper portion of the 
figure is an idealized sketch of the pressure time history showing chordwise 
propagation of a wave . (Note that in this sketch the leading and trailing 
edge positions are opposite to those of the memo scope time histories, but 
are in conformity with the contour plots . ) An event begins at the leading 
edge at t 0 and reaches a chordwise station X after some interval At (x) • The 
event repeats itself after one period, T, has elapsed. 

It was stated earlier that the constant altitude contour lines would 
form a nested family of curves, and this has been shown in previous figures. 
The locus of the. extrema of these contour lines will represent the direction 
of the ridge line on the X, t plane, and hence can be used to calculate the 
propagation rate. The lower portion of Fig. 13 shows a two-segment represen- 
tation of the locus of ridge line extremeties as it would be superimposed on 
one of the earlier contour plots. The point 0 is the extension of the locus 
to the leading edge and the point 2 is the extension to the trailing edge 
(recall that the memo scope photos cover the range 0.0119 — X — 0-90) . The 
point 1 represents the location of the kink in the locus, if it exists. 
Accordingly, the intervals At^ and At 2 represent the time required for the 
pressure wave to travel a distance AXp or a distance of one chord length, 
respectively. Two wave velocities, normalized with respect to the free- 
stream velocity, will be calculated from these parameters: the average 
normalized wave velocity for propagation over the entire chord, 

v = (7) 

W VAt 2 

and the initial normalized wave velocity for propagation over the leading 
edge ■ region up to the kink in the locus , 

c 

V = - 
W 1 V 

If no kink exists it is assumed that the initial and the average wave 
velocities are the same and only Eq. (j) is used. 

Some examples of the use of these procedures are found in Figs. Ik and 15 
for frequencies of 31 and 75 Hz , respectively. These are expanded memo scope 
plots in which less than a full period was displayed to yield greater accuracy 
in the region of dynamic stall cell propagation. The upper portions of these 
figures are again the time histories for 11, Ik, and l6 deg mean incidence 


AXl 


At 


1 „. 


(8) 
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angle, but the viewpoint is now different. Whereas in previous plots a 
three-dimensional effect was attained by viewing the pseudo-surface from 
some elevated position, the viewpoint in Figs. Ik and 15 is from ground 
level. In other words, only the leading edge pressure was referenced to the 
horizontal datum in previous figures while all time histories are referenced 
to the horizontal datum in Figs. lU and 15 • 

The lower portions of these figures are again the constant altitude 
contour plots, but these now have superposed ridge line loci as in the 
schematic of Fig. 13 • It is interesting to note that the time histories for 
f = 31 Hz in Fig. 14 have virtially no secondary pressure wave's and the 
corresponding ridge line loci are relatively straight, with only minor kinks. 
Conversely the time histories for f = 75 Hz in Fig. 15 all have both secondary 
and tertiary pressure waves that are quite prominent . For both a ^ = 11 deg 
and ofjyj = Ik deg the secondary wave seems to merge with the primary waVe and a 
pronounced kink is evident for the contour plots for these two cases. Note, 
however, that for = 16 deg the secondary wave remains distinct to the 
extent that many of the contour lines have two extrema and the locus of the 
primary ridge is not kinked. This lends credence to the assertion made 
earlier that the kinks appear to be associated with wave coalescence. 

Normalized wave velocities, both average and initial, were calculated 
for all frequencies over the range =11 deg toa M = 18 deg and are plotted 
in various forms in the next several figures. In Fig. 16 these wave velocities 
are plotted. versus airfoil reduced frequency and it is immediately evident 
that: l) all wave velocities are consistently less than one -half the free-stream 
velocity, 2) the initial wave velocities are consistently less than one-third 
the free-stream velocity, 3) the wave velocities at each mean angle of attack 
are reasonable well correlated (see next figure), and !<•) the wave velocities 
are strong functions of oscillatory blade frequency, varying almost directly 
with k. 

The correlation alluded to above is well depicted in Fig. 17 in which the 
wave velocities are cross-plotted versus mean incidence angle. With some 
exceptions these plots are generally horizontal, indicating that the wave 
propagation phenomenon is not a strong function of incidence angle, once the 
stall cell has broken away from the leading edge . Additional insight is 
afforded by defining two additional correlative parameters, the reduced 
wave velocity (average and initial), 

' *L W 1 ) 
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and the average and initial wave number 

•' * ' " k = l/y- 

• • w w . 

. " ’ k = l/v 

Z . . v w i 

With these definitions the wave propagation phenomenon can be well correlated 
within a relatively narrow band of parameter values as shown in Figs . 18 and 

19. ^ r; . ■ ' ; :r ; v: "- ' ^ 

Comparisons were made with results from other experimental programs 
(Refs. 8, ,20, 21, and 22) and the results are shown in Fig. 20. Pressure 
time histories presented in these references' were examined and the chordwise 
propagation,, of peak pressure was used to compute an approximate average 
chordwise. propagation rate for each case. A variety of airfoils were used 
in these, other studies, both thick and thin, including the NACA 0012, NACA 
0006, and two Boeing; prof iles . Other primary parameters also differed from 
those of the, UARL^ experiment, such as M, and a. All of these parameters 
are listed in the box at the top of the figure ,' together with the source of 
the data (reference and figure number), keyed to the symbol used in the plot. 
The present results from Fig. 1 6 (for the average normalized wave velocity) 

are indicated by, the envelope curves surrounding the grey area. 

if*- .£* '■* • rV. • ; ' .- - 

.• ... * ' * i , t * p f % ... 

r , Two impprtant results are contained in Fig. 20. The first is the 
obvious one r . that in all cases the wave velocity is considerably less than 
the,. free-stream,. velocity, and the frequency trend followed by most of the 
data frOT.other .sources is similar to that of' the present study. The second 
is revealed by noting that all data points lying below the envelope curve 
had amplitudes of approximately a = 5 deg, the data in the current study had 
amplitudes, of « =, 8' deg, and the single point from Ref. 20, which is above 
the envelope curve had an amplitude of a = 15 deg. It is true that these 
data also had significant Mach number and shape difference, but these 
differences, do . not, tend to separate the data as much as the amplitude 
differences^. Thus, it may be tentatively concluded that increasing the 
amplitude of motion tends to increase the wave velocity associated with 
separation/ ‘ . . .. .... ... . " „ .. ‘‘ 

Harmonic Data Analysis 

The first harmonic pressure amplitude, A-^(x), has been plotted versus X 
in Fig...21 for = 3> Uj lA, and 18 deg. Similar plots for Ag(x) are 
found in Fig. 22. For = 3 deg a single frequency, f = 50 Hz , is shown 
because of the lack of significant variation with f in both first and second 
harmonics (cf. Table II )• All frequencies are plotted for the other three 
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mean incidence angles. In Fig. 21 it is seen that an increase in generally 
causes a decrease in the first harmonic amplitude near the leading edge but - 
has little effect on A^ over the aft portion of the chord. In Fig . .22, the 
second harmonic amplitude increases from a negligibly small level at = 3 
deg to a measurable but still small value at higher incidence .angles .. ...The 
effect of increasing frequency is to increase the first harmonic amplitude 
for all cases examined (Fig. 21). A similar but less pronounced trend also 
occurs for the second harmonic amplitude over the aft portion of the airfoil, 
as shown in Fig* 22. However, an opposite trend in which an increase in f 
causes a decrease in second harmonic amplitude, occurs near the leading edge 
for intermediate values of Cross plotted data are shown in Figs . 23 and 

2k for the leading edge station for first and second harmonic pressure 
amplitude versus reduced frequency. Again it is seen that the first harmonic 
generally increases with k (Fig. 23) while the second harmonic generally 
decreases with k (Fig. 2k ) . The notable exceptions to the latter are = 3 
deg which is very small, and« M = 18 deg. 

Figures 25 and 2 6 contain representative chordwise variations in first 
and second harmonic pressure phase angles . Only the 50 Hz first harmonic 
data is shown for = 3 deg in Fig. 25 because other frequencies for this 
angle were shown previously in Fig. 6. It is seen that the primary ' effect ' r -' 

of increasing incidence angle on the first harmonic phase angle is to shift v 
the distribution from leading edge lag (negative phase angle ) at low or^ to 
leading edge lead (positive phase angle) at high ar^. In the latter case the 
phase angle distribution is strongly associated with the rearward propagation 
of the stall cell or pressure wave. The effect of increasing frequency at' ' 
low incidence, say = 11 deg, is to switch the behavior from a relatively 
level trend at f = 12.5 Hz to one with leading edge lag at f = 98.5 
which is a change towards potential flow. This trend is also seen at higher . 
values of or^, but here the chordwise phase angle distribution never reaches 
a condition of leading edge lag. 

In Fig.. 26 only the = 14 deg condition is shown for the second harmonic 
phase angle distribution. The second harmonic amplitudes are so small at low 
or jyj that the second harmonic phase angles are rendered meaningless, and the 
distributions at values of ar^ greater than Ik deg are much the -same as those _ 
at 14 deg. It is seen that there is a considerably wider phase variation 
over the chord than in the first harmonic distribution, and that the second 
harmonic wave at the leading edge leads the trailing edge, implying a rearward 
wave propagation here also. 
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Analysis of Ramp Cam Data 


A complete description of the ramp cam test is given in Ref . 17 , and only 
a brief discussion, pertinent to the analysis performed will be given here. 
Figure 27 shows the incidence angle schedule, relative to mean incidence, for 
the forward ramp cam. The upstroke region, CDAB, was designed to have a con- 
stant angular velocity (excluding the corner regions ) that is twice that of 
the downstroke region, BC. The cam was also operated backwards in which case 
the upstroke, region had an angular velocity that was one-half that of the 
downstroke region. This is summarized in Table V in which A = ca/2V, the 
angular velocity parameter, is tabulated for the cases studied in Ref. 17 • 

The nominal first harmonic frequency, f^, is the reciprocal of the period of 
one cycle of motion. (See Table I for a complete listing of the individual 
cases ruri.) 

It is seen in Table V that certain pairings of conditions have equal 
angular velocities. For example, a forward upstroke at f = 7 • 5 Hz and a 
backward upstroke at f = 1^.3 Hz both have A = 0.0020, and a forward down- 
stroke at f = 20 Hz and a backward downstroke at f = 10 Hz both have 
A = -0.0028. This design was specifically chosen to match angular velocities 
in these local regions during an upstroke (or a downstroke) while having 
large differences during a downstroke (or an upstroke) within one cycle of 
motion. This separate matching of upstroke and downstroke regions is shown 
schematically in. the upper and lower panels of Fig. 28 in which the solid line 
represents a forward motion and the dashed line represents a backward motion. 
(Note that these are idealizations of the actual motions which did not have 
sharp corners.) For upstroke matching, the forward motion frequency is half 
that of the backward motion and matching occurs over CDAB. For downstroke 
matching, the forward motion frequency is twice that of the backward motion 
and matching occurs over BC. In Ref. 17 these regions were examined using 
integrated normal force and moment loops as a basis of comparison. In 
general, it was found that equal angular velocities yielded nearly equal 
integrated results . 

A more detailed examination is possible using the pressure time histories, 
and for this purpose separate time domains were considered. First, the 
dimensionless upstroke time, t u , is defined over the upstroke region from 
the zero displacement point. A, to the maximum displacement point, B, 
normalized with respect to the time interval over AB (see Fig. 27). Similarly, 
the dimensionless downstroke time, tp, is defined over the downstroke region 
from maximum displacement, B, to minimum displacement, C, normalized with 
respect to the time interval over BC. In the latter case the entire downstroke 
region BC is considered because stalling, which usually occurs prior to the 
peak amplitude, continues over much of the downstroke region, followed by 
stall recovery near the end of the region. Conversely, only half the upstroke 
is needed because stall, does not usually occur until well beyond the mean 
incidence angle . 



Matching angular velocity regions fora^ = 11 and 14 deg are examined in 
detail for X = 0.0119 and O.2583 in Figs. 29 through 32. In each of these 
figures the solid lines are for forward motion and the dashed lines are for 
backward motion. In all cases the differences between forward and backward 
motions are greater for = 11 deg (Figs. 29 and 30 ) than for ajyj - lb deg 
(Figs. 31 and 32). In particular, Fig. 29 shows the results f or = 11 deg 
over the upstroke region AB of Fig. 28. For the lower angular velocity, 

A = 0.0020, the forward and backward pressure response is in reasonably good 
agreement, indicating the initiation of stall at the leading edge near 
t u = 0.5 to 0.6. Since there is a nearly linear angular variation this would 
correspond to an incidence angle of approximately 15 to l6 deg. A wider 
disagreement is found at A = 0.0028, with large differences in the time at 
which the various events occur for forward and backward motions. In Fig. 30 
for the downstroke region BC at = 11 deg there is a large discrepancy at 
A = -0.0020 but good agreement at A = -0.0028. Note that the flat region 
in the pressure at low t^ followed by a sudden rise and a gradual decline as 
tp increases represents the process of stall recovery (during the sharp rise) 
followed by decreasing pressure along a nearly linear characteristic as a 
decreases. Thus Fig. 30 indicates that at the lower (in magnitude) of the 
two angular velocities the backward motion experienced a less severe stall 
than did the forward motion. In contrast. Figs. 31 and 32 for the upstroke 
and downstroke regions at = 14 deg show very good agreement throughout . 

This would indicate that at = 11 deg the airfoil penetration into stall is 
not sufficiently deep to cause a large stall effect under all circumstances 
whereas at ay = 14 deg the additional 3 deg penetration into stall causes a 
more complete stall. It also appears that, for the limited range of parameters 
considered here, at low incidence the stalling behavior is dependent on past 
history (i.e., the large differences in angular velocity in the preceding 
upstroke or downstroke region) while at high incidence the stalling behavior 
is more dependent on local conditions such as the instantaneous value of A. 

Further clarification may be obtained from Figs . 33 and 3*+ which are the 
memo scope photographs for = 11 and 14 deg corresponding to eight of the 
cases just considered. In each figure the upper row is for forward motion 
while the bottom row is for backward motion. In the left columns the upstroke 
angular velocity has been matched and in the right columns the downstroke^ 
angular velocity has been matched (corresponding to the left hand panels in 
Figs. 29 through 32). In particular, consider the right column of Fig. 33 
which shows qualitatively the disagreement between the curves for A = -0.0020 
in Fig. 30* As stated earlier, the backward motion in this case does not 
stall so deeply that an abrupt recovery occurs while the forward motion does . 
Again, in contrast. Fig. 3^ shows that for = 14 deg there is good qualita- 
tive agreement between forward and backward motions in the appropriate matching 
regions. Indeed, there is generally similar behavior in all regions for this 
higher incidence angle, with a precipitous stall followed by a flat deep stall 
region, followed by an abrupt stall recovery. 

20 


CONCLUDING REMARKS 


As a result of this investigation additional insight has been gained in 
the study of the dynamic stall phenomenon. An ultimate solution and under- 
standing of this problem, however, will be obtained only when sufficient 
experimental data have been gathered and analyzed, and when more powerful 
theoretical techniques have been employed. Hence, the following list contains 
a set of tentative conclusions obtained from this study which are accurate 
and valid within the present framework, but which should be tested further as 
additional work is completed. First, for sinusoidal motion, the conclusions 
are : 

1. The dynamic stall phenomenon is characterized by a pressure wave 
which propagates from leading to trailing edge and is sometimes 
accompanied by secondary and tertiary waves. 

2. The average pressure wave velocity for rearward propagation over the 
entire chord is consistently less than one-half the free-stream 
velocity. 

3- The initial pressure wave velocity for rearward propagation in the 
vicinity of the leading edge is always less than or equal to the 
average wave velocity, and is consistently less than one-third the 
free-stream velocity. 

4. If the initial wave velocity is significantly different from the 

average wave velocity it is usually associated with the coalescence 
of a secondary pressure wave with the primary wave . 

5» The pressure wave velocity is a strong function of reduced frequency, 
varying directly with reduced frequency. 

6. The pressure wave velocity is a weak function of mean incidence 
angle and does not vary significantly with changes in incidence. 

7* Correlation with the results of other investigations indicates that 
increasing the amplitude of motion increases the wave velocity, but 
in all cases considered the wave velocity remains less than one-half 
the free-stream velocity. 

8. Near the leading edge the first harmonic pressure amplitude generally 
increases with increasing frequency while the second harmonic ampli- 
tude decreases . 



9* As mean incidence angle is increased the first harmonic chordwise 
pressure phase angle shifts from leading edge lag at low incidence 
to leading edge lead at high incidence . This trend is modified ty 
increasing frequency which causes the leading edge lead at high 
incidence to "be diminished. 

10 . At low incidence angles the results are in good agreement with 
potential theory. 

11 . At intermediate incidence angles, in a mixed flow condition, dynamic 
stall occurs near the end of the upstroke of the motion. Within the 
same cycle a potential flow condition exists over the beginning of 
the upstroke, which is terminated by stall inception. As the reduced 
frequency is increased the potential flow region persists to higher 
mean angles of attack. 

12 . An increase in frequency at high incidence generally decreases the 
severity of the dynamic stall, although stall recovery appears to 
be independent of frequency. 

For ramp cam motions the conclusions are: 

1. At high mean incidence angles in which stall penetration is deep, 
the stalling behavior is primarily dependent on local conditions such 
as the instantaneous value of the angular velocity of motion. 

2 . At intermediate mean incidence angles, in which stall penetration is 
less deep, the stalling behavior appears to be more dependent on 
past history than on local conditions . 


United Aircraft Research Laboratories 
United Aircraft Corporation 

East Hartford, Connecticut, September 28, 1973 
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TEST CONDITIONS ANALYZED 


































































TABLE II. - PRESSURE HARMONICS FOR SINUSOIDAL MOTION 
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- TABLE III. - CONCLUDED 
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TABLE V 


VALUES OF CONSTANT ANGULAR VELOCITY A 


FOR RAMP CAM TIME HISTORIES 
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-CYCLE SIGNAL-AVERAGED DATA 



DIMENSION LESS TIME, t* 

FIGURE 1. COMPARISON BETWEEN RAW DATA AND SIGNAL-AVERAGED DATA 





TIME, t, sec-. 


FIGURE 2. SAMPLE TIME HISTORY OF UNSTEADY PRESSURE DATA 





CHORD PRESSURE COEFFICIENT DIFFERENCE, ACp 


a M = 3 DEG 
f = 98.5 Hz 
k = 0.3308 



TIME, t, sec 

FIGURE 3. SCHEMATIC TO CLARIFY TIME HISTORY AIM'D CONTOUR PLOTS 
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DIMENSIONLESS CHORD, X 

FIGURE 4. COMPARISON OF HARMONIC AMPLITUDES FOR SIGNAL-AVERAGED 


DATA AND RAW DATA 








DIMENSIONLESS CHORD, X 


FIGURE 


COMPARISON OF HARMONIC PHASE ANGLES 
FOR SIGNAL-AVERAGED DATA AND RAW DATA 











alt) 
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FIGURE 9. PERSISTENCE OF POTENTIAL FLOW BEHAVIOR AT HIGH FREQUENCY 







FIGURE 10. EFFECT OF FREQUENCY ON PRESSURE TIME HISTORIES AND CONTOURS AT a M = 12 DEG 




s 
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FIGURE 11. EFFECT OF FREQUENCY ON PRESSURE TIME HISTORIES AND CONTOURS AT a m = 14 DEG 





ait) 



FIGURE 12. EFFECT OF FREQUENCY ON PRESSURE TIME HISTORIES AND CONTOURS AT a M = 18 DEG 
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FIGURE 15. TIME HISTORIES AND CONTOUR PLOTS WITH SUPERPOSED RIDGE LINE LOCI FOR f = 75 Hz 

























FIGURE 19. WAVE NUMBER AND REDUCED VELOCITY FOR INITIAL WAVE PROPAGATION 
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FIGURE 20. WAVE VELOCITY COMPARISONS 





























C PRESSURE AMP 


SYMBOL 
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HARMONIC PRESSURE AMPLITUDE, A 2 
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FIRST HARMONIC PHASE ANGLE, P,, DEG 













DIMENSIONLESS DOWNSTROKE TIME, t D 

1 

PERIOD, Ti = — 

f 1 


FIGURE 27. ANGULAR DISPLACEMENT TIME HISTORY FOR FORWARD RAMP MOTION 




FIGURE 28. RAMP MOTION SCHEMATIC SHOWING UPSTROKE AND DOWNSTROKE 
MATCHING REGIONS 










forward; fi>;i 4.20 Hz ;. — : forward^;; =19.74 Hz 
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FORWARD, f, =7.5S;Hz FORWARD, f, = 9.96 Hz 

BACKWARD, = 14.12 |Hz BACKWARD, f, = 19.73 Hz 
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FORWARD, f, = 14.19 Hz FORWARD, f , = 19.75 Hz 

BACKWARD, f, = 7.63 Hz BACKWARD, f, = 9.89 Hz 
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FORWARD, fj = 7.55jHz / ' FORWARD, f, = 14.19 ( Hz 
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FIGURE 34. COMPARATIVE PRESSURE TIME HISTORIES FOR FORWARD AND BACKWARD RAMP MOTION FOR a M = 14 DEG 
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